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Two parallel dielectric waveguides can exchange energy if the field 
carried by one guide reaches the other guide. We consider only the case 
of coupling between degenerate modes of dielectric waveguides (optical 
fibers, etc). Degenerate modes have equal phase velocities, but their 
transverse field distributions need not be identical. 

The coupling theory presented in this paper applies to dielectric 
waveguides of arbitrary shape and arbitrary distribution of refractive 
index. The dielectric media of the guides as well as the surrounding 
medium are allowed to be lossy. The coupling coefficient is obtained 
by means of perturbation theory. It is shown that whereas lossless 
degenerate modes can exchange their power completely, lossy modes 
tend to equalize their power. 

The theory is applied to the problem of crosstalk between cladded 
dielectric slab waveguides and cladded optical fibers embedded in a 
lossy medium. 

Since a lossy surrounding medium also causes an increase in the loss 
of the guided modes, formulas for this additional loss are presented. 
It is shown that additional mode loss results even for a lossless sur- 
rounding medium if n :i k > ft. (n :i = index of surrounding medium, k = 
free space propagation constant, ft = propagation constant of guided 
mode.) 

I. INTRODUCTION 

Optical fibers appear attractive as waveguides for the transmission 
of light. Since many such fibers arc likely to be bunched together to 
form a cable, the problem of crosstalk between the different fibers of 
the cable arises. There are several sources of crosstalk. Light scattered 
in one fiber may excite a guided mode of an adjacent fiber. This 
double scattering mechanism is discussed in the next paper in this issue. 
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The present paper is concerned with the direct coupling of power be- 
tween two parallel dielectric waveguides. Even though it addresses itself 
to crosstalk in optical fibers the theory is more general and applies to 
dielectric waveguides of arbitrary distribution of the refractive index. 
The modes carried by the two guides need not be identical. However, 
we assume that the phase velocities of the modes in either guide are 
identical. We call modes with identical phase velocity degenerate 
modes. We furthermore admit the possibility that the media of the 
waveguides as well as the surrounding medium in which the guides are 
embedded may be lossy. The restriction to degenerate modes is no 
serious limitation for cases of practical interest. Provided that the 
coupling coefficients are constant, appreciable coupling of power from 
a mode in one guide to a mode in the other guide is possible only if 
the modes have identical phase velocities. 

The coupling theory is based on finding the change of the propaga- 
tion constant of the mode of one guide caused by the presence of the 
other guide. This perturbation theory is directly based on Maxwell's 
equations and is thus independent of the composition and geometry 
of the coupled waveguides. 

We obtain the coupling coefficient in form of an integral over the 
cross section of the waveguides. The integral is taken over the scalar 
product of the electric field vectors of the two guides. The general 
coupling formula is then used to obtain the coupling between two 
cladded slab waveguides as well as the coupling between two cladded 
round optical fibers. The crosstalk between these two types of wave- 
guides is discussed. The crosstalk can be reduced by increasing the 
separation between the two guides and by increasing the loss of the 
medium in which the guides are embedded. Coupling between dielectric 
fibers has previously been treated by A. L. Jones 1 and by R. Van- 
clooster and P. Phariseau. 2 However, these authors consider only 
uncladded lossless fibers. 

II. COUPLED LINE EQUATIONS 

The general properties of coupled waveguides can be studied with 
the help of coupled line equations. Neglecting the possibility of 
coupling to modes traveling in the opposite direction we can write the 
general coupled line equations of two degenerate modes as follows: 3 

4£= -&A+ Cl B, (la) 

az 
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— = -tfB + CiA. (lb) 

az 

The propagation constant (3 is assumed to be the same for either 
mode; A is the amplitude of the mode of one guide and B that of the 
mode of the other guide. If each guide can support more than one mode 
these modes would have to be included in the coupled line equations. 
However, it is well known that only modes with equal phase velocity 
exchange a significant amount of energy if they are coupled together 
by a length-independent coupling mechanism. The restriction to only 
two modes is thus an excellent approximation. The coupling coefficients 
Ci and c-> are allowed to be different since the two modes can be 
different even though they are degenerate. However, if the waveguides 
are lossless, conservation of power imposes the condition 8 

c, = -c*. (2) 

The asterisk indicates complex conjugation. 

In the following analysis we assume that /? as well as C\ and &> may 
be complex quantities since we want to include the case of lossy modes. 
Equation (2) is thus not assumed to hold exactly. Since /? as well as 
Ci and c-> are constants we can immediately solve the coupled line 
equations and obtain solutions in the form 

A{z) = \ {A (e-"" + •"*) + (l)* B (e-"* - «"*)}»-*■, (3a) 

B(z) = \ {b (<T , ' 4 * + *"*) + (£f A„(e- A * - e "*)}<f '* (3b) 

with 

A/3 = iVe& . (4) 

The coefficients A anil B„ are the field amplitudes A and B at z = 0. 

It is apparent that equations (3a I and (3b) are composed of the 
superposition of two new normal modes with propagation constants 

A = + A£ (5) 

and 

2 = & - A/3. (6) 

For complex values of C\ and c- we obtain also a complex value of the 
change of the propagation constant A/?. In that case, it is clear that 
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the two new normal modes do not only have slightly different propaga- 
tion velocities but also suffer different amounts of loss. 

For the study of crosstalk we need the solution with B n = 0. That 
means we assume that at z = only mode A is excited. If we restrict 
ourselves to the case 

| A/3z | « 1 (?) 

we find that the ratio of the amplitude B to amplitude A for distances 
z = L for which (7) is valid is given by 

The absolute square of the ratio of B/A is the power ratio at z = L 
in the two guides. Using the letter V for this power ratio we obtain the 
following useful formula for the distance L at which the power ratio 
in the two guides is V : 

(y)! (9) 



L = 



I A/3 | 

For low-loss modes equation (2) must hold at least approximately so 
that we obtain from (9) in the low-loss case 

L ' tSt- (10) 

The low-loss assumption refers only to the actual loss suffered by the 
modes propagating inside of the guides. The loss of the surrounding 
medium in which the guides are embedded can be arbitrarily high. 

Finally we consider two special cases. If we assume that the two 
guides are lossless and assume again B„ = we obtain from equa- 
tion (3) 

A(z) = A coB(Afc)e-"", (Ha) 

B(z) = -i(^jA sm (A/te)e""". (lib) 

In the lossless case assumed here A/3 is real and it is clear that the 
two guides exchange power after a distance D given by 

D - w (12) 

If, on the other hand, A/3 is complex, as it would be in the case of lossy 
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guides, and if we make z .sufficiently large that the exponential factor 
cxp(iAfiz) is much less than unity, we obtain again for B„ = from 
equation (3) 



A® = \A+ mm (13a) 



' iCa • -i(/J+A0)i 



and 



B ® = %<7 A « ■ ( 13b ) 

Since the absolute value of the ratio c- 2 /ci is nearly unity we see that 
the power carried by the two guides equalizes in the lossy case pro- 
vided that we let both modes travel far enough. This power equaliza- 
tion may occur after many power exchange cycles or it may occur 
even before one cycle is completed depending on the magnitude of 
the imaginary part of A/?. Of course both modes suffer additional loss 
because of the complex value of the propagation constant /3. The 
magnitude of /3 is much larger than that of A/3 for this discussion to 
be meaningful. 

The preceding discussion shows that the crosstalk between the two 
coupled guides is determined by the change in propagation constant 
A/3 and not by the individual coupling constants Ci and c 2 . It is thus 
sufficient to determine A/3. 

III. MODE COUPLING THEORY 

We assume that two dielectric waveguides are positioned parallel 
to each other. Each guide is a cylindrical structure with a refractive 
index distribution that is independent of the z coordinate (z is the 
direction parallel to the axis of the structure) but which can have an 
arbitrary dependence on the transverse x and y coordinates. However, 
it is assumed that the index distribution is such that it produces a 
dielectric waveguide capable of supporting guided modes. We assign 
a refractive index distribution Hj to guide one assuming that ni has 
the constant value of the refractive index ?2 : , of the background ma- 
terial in the region occupied by guide two. A similar assumption 
applies to the index distribution of guide two. It is assumed to have 
the constant value n& of the background in the region occupied by 
guide one. The square of the index distributions ?u and /?-_> is shown in 
Fig. 1. Since it is the square of the refractive index that enters 
Maxwell's equations we use the following expression 

n 2 = (n 2 - nt) + (n\ - n\) + n\ . (14) 
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Fig. 1— Distribution of the square of the refractive index used to describe the 
two dielectric waveguides. 

Owing to the definitions of n x {x, y) and n 2 (x, y) it is clear that (14) 
becomes n = n\ in the region of guide one, n = n\ in the region of 
guide two, and n = n\ outside the regions of either guide. The re- 
fractive indices rii , n 2 , and n 3 need not be real quantities but are 
allowed to be complex if the media are lossy. 
The starting point of the perturbation theory is Maxwell's equations 



and 



V X H = iue n 2 E 



V X E = -iwMoH. 



(15a) 



(15b) 



The time dependence of the fields was assumed to be given as exp (iut). 
The constants e and mo are the permittivity and permeability of the 
vacuum and the square of the refractive index is given by (14). Next 
we introduce the perturbation assumption that the electric field E 
and the magnetic field H are very nearly given by the superposition 
of the mode fields of each guide in the absence of the other. We thus 
write 



and 



E - dBi + bE 2 + e 



H = oH, + bH 2 + h. 



(16) 



(17) 



The coefficients a and b are being determined by the theory. The 
field E„ , H v (i/ = 1, 2) being a mode of guide v has the z dependence 
exp (— 1*/3 2). It is assumed that the propagation constant is the same 
for the modes of either one of the guides in the absence of the other. 
As indicated in Section I, such modes are here called degenerate. The 



COUPLING OF DEGENERATE MODES 1797 

electric and magnetic fields of the unperturbed modes of either guide 
satisfy the equations 

V, X H„ - #o(2 X H„) - iue a n 2 ,E, = v - 1, 2 (18a) 

and 

V, X E, - i/3 (z X E„) + icoMoH, = 0. (18b) 

z is a unit vector in z direction. The symbol V, indicates the transverse 
part of the V operator. The electric vector e and the magnetic vector h 
appearing in (16) and (17) are assumed to be small perturbations 
that are needed to express the change in the field shape resulting from 
the close proximity of the two guides. In addition to the changes in 
the field shape we must allow for the possibility of a change in the 
propagation constant which can be assumed to be small if the coupling 
of the two guides is sufficiently loose. We thus assume that the new 
fields E and H of (16) and (17) have a z dependence of the form 
exp (—i&z) with 

|8 = /3 + A/3. (19) 

We have seen in the last section that the change A/3 of the propagation 
constant determines the coupling of the two guides. It is thus the 
quantity we are most interested in. 

We substitute (16) and (17) into (15a) and (15b) using (14). After 
cancellation of a number of terms with the help of (18a) and (18b) we 
are left with the following system of equations: 

V, X h - ifi (£ X h) - iu>e [{n\ - n\) + {n\ - n\) + n\]e 

= iA0(az X H, + bz X H 2 ) 

+ iu>e [(nl - ?i 2 3 )oE, + (n] - r&)bE 2 ] (20) 

and 

V, X e - ij8 @ X e) + iw Mo h = /A/3(a| X E, + bz X E 2 ). (21) 

Products of A/3 with e and h have been omitted from (20) and (21) 
since they are small of second order. 

Since we are interested in obtaining A/3 from our theory we must 
eliminate the dependence on the spatial coordinates from (20) and 
(21). Concentrating for the moment on the left-hand sides of these 
equations we obtain by scalar multiplication of (20) with E, and of 
(21) with H! and by subtracting the two equations and integration 
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over the entire cross-sectional area 

f {Ex-tV. X h - i(3 (z X h) - «we ((n a , - n 3 ) + (n\ - n\) + n 2 3 )e] 

- Hi-[Vi X e - i/3„(z X e) + w/ioh]] dte cfy. (22) 

The terms inside of the square brackets are small of first order. We 
now assume that the field intensity of Ex in the region of guide two 
has decayed sufficiently to render it also small of first order. Since 
n\ - n\ is different from zero only in a small area in the vicinity of 
guide two the product of E, with (>?* - n 2 3 )e is a term that is small 
of second order. In first order of perturbation theory this term can 
be neglected. Performing a partial integration on the terms containing 
the V, operator and rearranging the terms of the mixed products 
allows us to write 

j {h-[V, X E, + iM X BO - ^MoHJ 

- e-[V ( X H, + iM X HO + iWflE,] } dx dy. (23) 

Comparison of (23) with (18) shows that this expression would vanish 
if instead of Ei and H x we had used the complementary fields E~ and 
H"[ that result from the original field if ft, and w are replaced by -ft, 
and —w. We have thus shown that by scalar multiplication with E~[ 
and H~, subtraction, and integration over the entire cross-sectional 
area the field terms e and h can be made to vanish from (20) and (21). 
A similar procedure can, of course, be carried through using E ~ a and H ~ a . 
We thus obtain from (20) and (21) the following set of equations: 



a A^-|(H7 



X E, + H x X ET) dx dy 

+ ue f (n 2 2 - n^ET-E, dx dy j 



+ 6 A&-J (H.\ X E 2 + H 2 X E7) dx dy 

+ a.€ f {n\ - nl)E- -E 2 dx dyj = (24a) 



and 



a A&-J (Ha X E, + H, X E" 2 ) dx dy 
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+ toe J {n\ - nI)Ea-Ej dx dy\ 

+ b\ m- J (Ha X E 2 + H 2 X Ea) dx dy 

+ coe J (»? - nJDE;-B a dx dy =0. (24b) 

Equations (24a) and (24b) still contain a number of second-order 
terms. Products of field terms with different indices are small of first 
order since the modes of the two guides are supposed to overlap only 
slightly. Products of field terms with different indices that are mul- 
tiplied by the first-order quantity A/3 can thus be neglected as being 
of second order. The term (n\ — n£)Ei -Ei is also of second order since 
it involves the square of the amplitude of Ei at the location of the 
opposite guide. Neglecting it and a similar term results in a consistent 
equation system containing only first-order terms. 

aA&-f (H7 X E, + H, X EI) dx dy 

+ bue I (n\ — ri3)E"i-E 2 dx dy = (25a) 
and 

fl^co / (^2 — n^Ea-Ei dx dy 

+ bA0z- f (Ha X E 2 + H 2 X Ea) dx dy = 0. (25b) 

The equation system (25) allows us to determine the coefficients a 
and b. Since we have a homogeneous system of equations we must 
require that the system detenninant vanishes. This latter condition 
leads to a determination of A/3. 

A/3 = ±co€ 

J (nl - nl)E~-E 2 dx dy J {n\ - n&Bi-B, dx dy 



j z- (E, X H7+E7 X H.) dx dy j i- (E 2 X Ha+E^ X H 2 ) dx dy 



(26) 
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The ratio of a/b is of no interest to us so that we do not need to write 
down the solution of (25) . 

Equation (26) is the final result of our coupling theory. In the general 
case of lossy dielectric media A/3 is a complex quantity. Equation (26) 
is very general and holds for any two dielectric waveguides carrying 
degenerate modes. The field expressions E~; and H~ are obtained by 
changing the sign of /3 and o> in the regular expressions for E, and H,- . 
In the special case of lossless media E~ and H^ become simply the 
complex conjugates of the original fields. 

'I for real n, . (27) 

H7 = HHJ 

Equation (26) can be simplified in another special case. If both wave- 
guides are identical, carrying the same mode, the two types of integrals 
appearing in the numerator as well as in the denominator of (26) 
become identical for reasons of symmetry. The expression for identical 
modes of two coupled identical waveguides can thus be written more 
simply 

f (nl - wDEVE, dxdy 

A/3 = ±a* — • (28) 

J |.(E, X H7 + E; X HO dx dy 

For real refractive indices the denominator of (28) is equal to AP, P 
being the power carried by the mode of one of the guides. It is neces- 
sary to normalize the mode fields of each guide so that field 1 as well 
as field 2 carry equal power P. This normalization was already 
implied in converting (26) to (28). Even if the medium surrounding 
the two waveguides is lossy it is possible to use the relations (27) and 
the identification of the denominator of (28) as 4P provided that the 
surrounding medium has a very small influence on the mode fields 
carried by each guide. It is thus permissible to use (28) in the form 

A/3 = ±|p / (nl - nDE* 2 -E, dx dy (29) 

in the special case of lossless dielectric media or, at least approxi- 
mately, even in the case that the surrounding medium is lossy if 
only the fields are so well guided that only a very small amount of 
power reaches the lossy surrounding medium. 
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IV. APPLICATION TO CLADDED SLAB WAVEGUIDE 

The dielectric slab is the simplest type of dielectric waveguide. 
Because of its simplicity much insight into the performance of dielec- 
tric waveguides can be gained by studying the slab waveguide. We 
consider two identical cladded slab waveguides as shown in Fig. 2a. 
Instead of solving the mode problem of the cladded slab waveguide 
exactly, we assume that the cladding is thick enough so that only a 
small amount of power reaches the surrounding medium. This assump- 
tion allows us to obtain approximate field solutions that are much 
simpler than the exact solutions. This case is, moreover, of the most 
practical interest since the cladding is meant to protect the guided 
mode from the disturbing influence of the surroundings. Our approxi- 
mation thus coincides with the case of real practical interest. Without 
going into the details of the calculation we state only the results. The 
field of one of the two guides can be expressed as follows: (For 
simplicity the field of one guide is expressed in the coordinate system 
shown in Fig. 2b.) 



E„ = 



A cos kx 
Be 11 + Ce~ yz 
[Fe~ pz 



^ x £ d 

d ^ x ^ D 
D < x < oo . 



(30) 



The field is continued symmetrically for x < 0. Equation (30) gives 
the E component of the symmetric TE modes of the slab waveguide. 
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Fig. 2a — Sketch of two coupled cladded slab waveguides. 
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Fig. 2b— Sketch of a single slab waveguide used for equation (30). 
The other field components are obtained from E„ by the relations 



and 






O)/i OX 



(3D 



(32) 



The dependence of the field on time t and on the z coordinate is, 
as always, given by the factor 



.Hvt-Bt) 



e— - (33) 

which is being omitted from the equations. The parameters appearing 
in (30) are defined as follows: (k = co\/«oMo) 

« = (n 2 c k 2 - /3 2 )*, (34) 

y = (/3 2 - n 2 m k 2 )\ (35) 

P = (0 2 - n\k 2 )K (36) 

n c , n m , and n> A are the refractive indices of core, cladding, and outer 

medium. The boundary conditions lead to relations between the 

constants A, B, C, and F: 



r _ A K 7 ~ P f -T< r -aTCP-» 

2 (k + 7 ) 7 + P 

C = Ae yd cos Kd, 



(37) 
(38) 
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*" - ( X ^rr x je> D e-^ D - d \ (39) 

(7 + p)(k + 7 ) 
The coefficient A can be related to the power P carried by the mode, 

2o> M P 



A = 



Pd + V- 
y. 



(40) 



The coefficient F depends exponentially on the factor —y(D — d). For 
increasing cladding thickness this factor becomes very small. However, 
the factor B vanishes even more rapidly since it depends on the square 
of the same exponential decay factor. This shows that the field inside 
of core and cladding is altered only very slightly if the cladding is 
sufficiently thick. The eigenvalue equation differs from the eigenvalue 
equation for infinite cladding only by a term that also vanishes as 
the square of the aforementioned exponential decay factor. To a good 
approximation we are thus justified to use the eigenvalue equation 
for the infinite cladding guide. 

tan Kd = X - (41) 

K 

Once the mode fields are known it is an easy matter to evaluate 
the coupling coefficient (29). We obtain 

A 2 3 

(R is the distance between the core centers of the two slab wave- 
guides.) In case that the cladding material is identical to the sur- 
rounding medium we have p = y and (42) reduces to the expression for 
coupled slab waveguides to be found in N. S. Kapany's book. 4 

The same mechanism that is responsible for the coupling between 
the two slab waveguides also causes loss to the modes in each of the 
guides. Coupling occurs because some of the field in one guide reaches 
the region of the other guide. However, fields reaching out into the 
lossy medium between the guides also cause power loss to the mode 
traveling inside of the guide. Only if the surrounding medium is truly 
lossless can there be no additional loss to the guided mode. I calcu- 
lated the mode loss contributed by the surrounding medium assuming 
that it is infinitely extended and also ignoring the presence of the 
other guide. The power loss can be obtained to a good approximation 
by computing the power flow in transverse direction to the guide 
axis at the boundary between the cladding and the surrounding 
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medium. The validity of this approach hinges again on the fact that 
most of the power is carried inside of core and cladding. The approxi- 
mate calculation leads to the following expression for the power loss 
2« of the mode (a is the amplitude loss coefficient) . 



, _ 8«V Im (p)<T 2 " 

(3(l+yd)(K 2 + 7 2 ) |7+ P 



2a = , ** 7 *? W * ■ ; n- (43) 



It is assumed that n c and n m are both real; ?i :i is allowed to assume 
complex values. This loss expression has two interesting features. It 
shows that the loss decreases exponentially with increasing cladding 
thickness. It also shows that the loss is dependent only on the imagi- 
nary part Im(p) of p. For p to be real, two conditions must be satis- 
fied. First of all, the refractive index of the surrounding medium, n 3 , 
must be purely real. This means that the mode suffers loss whenever 
the surrounding medium is lossy. However, even if the surrounding 
medium is lossless, it is still possible for P to be purely imaginary. 
This happens when n :i k > (3. In this case there is no total internal 
reflection between the cladding and the surrounding medium and the 
evanescent field tail reaching into the surrounding medium does no 
longer decay exponentially but begins to radiate. Power is thus lost 
from the guided mode even though all the materials of the guide and 
the surrounding medium are lossless. 

Finally, we state the result of the coupling calculation and mode 
loss for the cladded slab waveguide carrying the symmetric TM mode. 
For the details of the field distribution of this mode in an infinite 
cladding guide see Ref. 5. The result of the evaluation of (29) in this 
case is 

AO n 2 c nyy 2 [2e {y ' pUD ~ d) " lfc-^-V"-" 1 " , , 

A/3 " " $[«* + nh*hd + ntf + 7 2 )] l } 

The power loss caused by the lossy surrounding medium is 

O 2 * I ,*, I 4 2 3 t„ ( P \-2y (D-d) 

8n c n m I n 3 I k 7 lm I -2 Je 

2a ~ Kyd(nU 2 + nh 2 ) + n c n 2 m (K 2 + 7 2 )] I n 2 mP + nh T K } 
It is apparent that the coupling expression (42) for the TE mode 
differs from that for the TM mode. The coupling as well as the loss, 
thus turn out to be polarization dependent. 

V. COUPLING OF THE HE n MODES OF ROUND CLADDED FIBERS 

In complete analogy to the cladded slab waveguide we now proceed 
to a discussion of coupling between HE n modes in two round cladded 
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Fig. 3 — Cross section of two coupled cladded fibers. 

fibers, Fig. 3. Again, we assume that the cladding is sufficiently thick 
to allow us to treat the field inside of core and cladding as being 
affected only very little by the presence of the surrounding medium. 
With this approximation in mind, we state the relevant equations for 
the mode field of a cladded fiber. 

The field inside of the fiber core is described by the set of equa- 
tions 



^ r ^ a, 



(45) 



a ^ r^ b, (46) 



E c = AJ x (nr) cos<p 
H, = BJ x (nr) sin tf> 
the field inside of the cladding is given by 

E z = [C//, (1 W) + DH[ 2 \iyr)] cos0 
H z = [FH^Xiyr) + GH\ 2 \iyr)] sin </> 
and the field in the surrounding medium is 

E t = MH[ l \i P r) cos^l 
H, = NH^iipr) sin ) 

The parameters k, y and p are given by the equations (34), (35), and 
(36) . The other field components are obtained from the longitudinal 
components by differentiation 

n i (. dE z . 1 dH,\ 



b £ r < » . 



(47) 






co/x 



dr 



(48) 
(49) 
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H r = - 



H A — — 






we 



+ coe n 



dE, 
d<f> 

2 dE_ 
dr 



(50) 



(51) 



The symbol r stands either for k, iy, or ip depending on the medium 
in which the field is evaluated. The same is true for n, it represents 
either n c , n m , or % . The boundary conditions establish the following 
relations between the amplitude coefficients: 



B 
A 



kan 



W + y) 



u ri 

U \-i<a 



J (ko) "[ 



+P 



\_ya H x (iya) J J 



N = 



c 


F 


Ji(ko) 


A 


B 


~ Hl l) (fya) 




D« 


iGfaQ, 



(53) 
(54) 



17 = ^ I . 



(p-y)b 



• c 



ff(i-$b-i)-*t*+x<+*i 



+ 1 (nl - nl) 
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The symbols «/ and «/i indicate Bessel functions of zero and first order; 
H ll) (iya) and ffi (1) (iya) are Hankel functions of zero and first order 
and of the first kind. We use here the notation of a Hankel function 
with imaginary argument used by Jahnke and Emde. 7 The coefficients 
M and N given by (55) and (56) do not contain any Hankel functions 
because the approximation for large arguments was used. The relation 
between the power P carried by the HEn mode in the waveguide and 
the amplitude coefficient A is given by the equation 
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P = I {^ t(a K ) 2 (J^(«a) + J 2 M) - 2J\{Ka){n\ + * p) 

H.^lfA^.t^q^^, (57) 

For sufficiently large cladding radius 6 the propagation constant of 
the mode can still be obtained from the eigenvalue equation for a rod 
embedded in a medium with refractive index n m . 

or 2 feM 1 ^ + J iH?\iya) _ t \\li a£ (j n ( K a) 1 



k \t/ 1 (ica) Kd/ m \ Hi (iya) I)\ m K XJ^kcl) ica 






(n'-n 2 )^]*- (58) 



This collection of formulas is sufficient to calculate the coupling term 
(29) for the coupling of the HE n modes in two parallel cladded 
optical fibers. We distinguish two cases. We assume that the modes 
are either both polarized perpendicular to the plane connecting the 
two fiber core centers or that they are polarized parallel to this plane. 
The coupling is different in the two cases, corresponding to the dif- 
ference in coupling between the TE modes and the TM modes of the 
slab waveguide treated in the preceding section. It must be pointed out 
that a perpendicularly polarized mode does not couple with a horizon- 
tally polarized mode. The two-mode assumption is thus still valid. 

The integral occurring in (29) cannot be solved generally for two 
coupled HEn modes. However, if we assume that the distance R 
separating the two cores is large compared to the core radius "a," we 
are justified in expanding the expression for the radial coordinate r 
(reaching from the center of guide one to the vicinity of guide two) 
in the following way (Fig. 3) , 



r, = VR 2 + r 2 2 + 2r 2 R cos 6 « R + r 2 cos 6. (59) 

The angle 8 describes the departure of the direction of r^ from the 
line connecting the core centers, r-> is the distance from the center of 
guide two to the endpoint of r\ . When the expansion (59) is used, 
the angular integration can be performed leading to Bessel functions 
of the imaginary argument ipr-> . The remaining integration over r 2 
now involves only products of cylinder functions and can be carried 
out. We obtain the following result for the case of horizontal polariza- 
tion: 
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For the case of vertically polarized modes we get 



(60) 
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Finally, we present also the expression for the power loss 2« of the 
HEn mode caused by the lossy surrounding medium. The calculation 
follows the idea described in connection with the slab waveguide. 
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(62) 
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The coefficients M and N are given by (55) and (56) , p r is the real 
part of p, and the symbol Im ( ) indicates the imaginary part of the 
expression in parenthesis (or brackets). 



VI. NUMERICAL RESULTS 

To illustrate the consequences of the theory presented in this paper 
a few sample cases have been computed. For comparison purposes 
I have checked the results of my theory against the curves for the 
coupling of uncladded fibers published by Jones. 1 The agreement 
between his theory and similar curves obtained from (60) and (61) is 
excellent. 

The amount of coupling between two optical fibers depends on the 
type of mode they are carrying, on the separation between the fiber 
cores, and on the loss of the surrounding medium in which the fibers 
are embedded. Figure 4 shows a set of curves describing the coupling 
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parameter a |A/3| as a function of the loss in the surrounding medium. 
The abscissa labeled "loss" measured in db indicates the plane wave 
loss in going through the surrounding medium for a distance R — 2b; 
that means it is the loss that is encountered in going perpendicular 
to the fiber axis from the cladding boundary of one fiber to the other. 
There are two curve parameters in Fig. 4. The two sets of curves differ 
by the cladding thickness normalized with respect to the core radius, 
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Fig. 4 — Coupling of round fibers. a|A/J| is the coupling parameter of equations 
(60) and (61), the abscissa indicates the loss (in db) that a plane wave would 
experience in going from the boundary of one fiber to the boundary of the other 
fiber. The solid lines hold for horizontal polarization, equation (60), while the 
dotted lines apply to the case of vertical polarization, equation (61). n m ka = 16, 
ya = 1.594, n c /n m — 1.01, Re(no) = n» . 
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(b - a) /a. The coupling is, of course, reduced for increasing cladding 
thickness since the fiber cores are farther removed from each other. 
The second parameter is the separation between the fiber claddings 
normalized with respect to the fiber core radius, (R - 2b) /a. Increas- 
ing values of this parameter again cause the fiber cores to be sepa- 
rated farther from each other. The loss values that form the abscissa 
hold regardless of fiber separation, always indicating the total loss 
that exists between the fiber claddings on their closest approach. The 
curves show clearly the decrease in coupling with increasing loss in 
the surrounding medium. The solid curves apply to the horizontal 
polarization of the two modes and are obtained from (60) . The dotted 
curves apply to vertical polarization and were computed from (61). 
Figure 4 holds, of course, only for one particular pair of round optical 
fibers. It was assumed that n m ka = 16, and that the ratio of core index 
to cladding index is n r /n m = 1.01. For simplicity, it was furthermore 
assumed that the real part of the index n :i is identical to the index n*, 
of the cladding material. The imaginary part of n 3 is varied in order 
to obtain the variable loss values of the abscissa. The abscissa can 
be expressed by 8.68 Im (n 3 )fc(J2 - 26). 

With the help of ( 10) , Fig. 4 can be used to evaluate the crosstalk 
between the two optical fibers to which the conditions of Fig. 4 
apply. Let us assume that we can tolerate a power ratio V = 10 -6 at 
the end, z - L, of the two fibers. If we use, for example, a = 1 /mi and 
allow L = 1 km we find that we must not exceed 

a | A/3 | ^ 10" 12 . (63) 

Any value appearing in Fig. 4 below this line thus leads to acceptable 
crosstalk levels. It is apparent that all curves with (6 - a) /a = 10 
are acceptable. For (6 — a) /a = 5 a certain amount of loss is re- 
quired to reduce the coupling and consequently the crosstalk below the 
desired level. 

Figure 5 presents the same data as Fig. 4 applied to two identical 
slab waveguides. It was assumed that the refractive indices as well 
as the guide dimensions of the slab waveguides correspond exactly 
to the round optical fibers of Fig. 4. (a = d, b = D). However, the 
two slab waveguides are less tightly coupled than the corresponding 
round fibers. The reason for this behaviour is the difference in the 
decay parameter ya (or yd) . For the round fibers, we obtained from 
the eigenvalue equation (58) the value y a = 1.594 while (41) results 
in yd = 1.997 for the slab waveguides. A larger value of the decay 
parameter indicates that exponential decrease of the field amplitude 
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outside of the fiber core is more rapid. The more rapid decrease of the 
field of the slab waveguide leads to less coupling. The solid and dotted 
lines have the same meaning as in Fig. 4. The solid lines apply to the 
TM mode which is polarized horizontally while the dotted lines apply 
to the vertically polarized TE mode. For simplicity we have assumed 
that the values of yd are identical for the two modes which, strictly 
speaking, is not quite true. However, we see from Fig. 4, as well as 
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Fig. 5 — Coupling of slab waveguides. Coupling parameter and loss have the 
same meaning as in Fig. 4. The solid lines indicate the TM mode (horizontal 
polarization) while the doited lines apply to the TE mode (vertical polarization). 
timka = 16, 70 = 1-997, n 8 /n. = 1.01. 
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from Fig. 5, that the horizontally polarized modes are coupled more 
tightly than the vertically polarized modes. 

Figure 6 presents a comparison between the coupling strength of 
round fibers and slab waveguides. Since equal geometry does not lead 
to identical values of the decay parameter we have arbitrarily used 
equal values of y a = yd = 1.594 for both types of waveguide. The 
solid lines indicate the results for the round fiber while the dotted 
lines apply to the slab waveguides. It is apparent that the slab wave- 
guides are more tightly coupled when the conditions are such that 
both guides have equal decay parameters. This appears reasonable if 
we consider that the slab waveguides have constant distance from 
each other over an infinite length while the round fibers have a point 
of closest approach so that the total amount of field overlap is less 
in this case. 

The lossy medium does not only serve to reduce the coupling be- 
tween the waveguides but it also has the adverse effect of increasing 
the loss of the modes traveling in the guides. Figure 7 is a plot of 
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Fig. 6— Comparison between the coupling of round fibers and slab waveguides 
(a = d, b = D) . The parameter value ya = 1.594 is correct for the round fibers. 
It has also been used for the slab waveguides to enable a realistic comparison. 
n m ka = 16, nc/n m = 1.01. 
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mode power loss versus plane wave power loss 2/c[Im (n : ,)a] in the 
surrounding medium obtained from (43), (44b), and (62). Both types 
of waveguides are represented in this figure. The loss values of abscissa 
and ordinate are both expressed in db. The decay parameters ya and 
yd (we use a — d and b = D in Fig. 7) are solutions of the respective 
eigenvalue equations and arc not adjusted arbitrarily. We have seen 
that the coupling of the slab waveguides was weaker than the coupling 
of the fibers under these conditions. It is thus not surprising to see 
from Fig. 7 that the losses of the fiber are larger than the slab wave- 
guide losses. It is interesting to see that the symmetric TM mode of 
the slab waveguide is slightly lossier than the symmetric TE mode. 

Assuming again that a = d = 1 /tm and limiting the allowable loss 
caused by the lossy surrounding medium to 1 db/km, we find that we 
can accept those curves with losses of less than 2aa = 10~ n db. Since 
the loss curves are nearly independent of the actual value of the loss 
in the surrounding medium (after an initial rise from zero mode loss in 
a lossless surrounding medium) we can exclude all curves above the 
10~ n line as unacceptable in case of a lossy surrounding medium. We 
thus see that the slab waveguide works well under all conditions shown. 
For the more important round fiber we must exclude the case of a 
cladding with (b — a) /a = 5. This means that the curves of Fig. 4 that 
lead to excessive coupling in the absence of sufficient loss in the sur- 
rounding medium are also unacceptable from the point of view of 
additional mode loss. We see that additional mode loss and mode 
coupling have similar trends. When the coupling between the fibers 
is too large we also have to contend with additional mode loss in 
excess of what we want to tolerate. A guide with low mode loss on 
the other hand is also sufficiently protected by its cladding so that 
coupling between adjacent fibers is not critical. 

The two remaining curves, Figs. 8 and 9, show the coupling param- 
eter and the additional mode loss for the case of a guide with lower 
core-to-cladding ratio. Only the case of the higher coupling for hori- 
zontally polarized modes has been plotted in Fig. 8. The remarks about 
coupling strength and additional mode loss hold true also for these 
conditions as can be seen from the figures. 

VII. CONCLUSIONS 

We have presented a perturbation theory for the coupling between 
two arbitrary dielectric waveguides made of lossy materials and 
embedded in a lossy environment. Only the case of two degenerate 
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Fig. 7— Additional mode loss caused by the lossy surrounding medium. The 
abscissa indicates the loss (in db) that a plane wave would suffer by traveling 
a distance equal to the core radius in the surrounding medium. n m ka = 16, 
n e /n m = 1.01. 

modes has been treated. The general theory has been applied to cladded 
slab waveguides and cladded round optical waveguides (fibers) with 
lossless cladding material but in the presence of lossy surroundings. 
We have found that the coupling between fibers depends very strongly 
on the separation between the two fiber cores. Far separated cores 
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result in loose coupling. The coupling between the fibers can be re- 
duced by increasing the losses of the surrounding medium. However, 
a lossy surrounding medium introduces losses to the waveguide modes 
to such an extent that the mode losses are unacceptably high if the 
loss in the surrounding medium is necessary to uncouple the fibers. 
From this point of view it appears as though losses in the surrounding 
medium are unsuitable to achieve uncoupling or reduction of cross- 
talk between optical fibers. However, the coupling mechanism dis- 
cussed in this paper is not the only contributor to crosstalk. Imperfec- 
tions in the fibers cause light scattering that can also be a source of 
crosstalk. 8 This mechanism is less strongly dependent on the separa- 
tion of the fiber cores so that it may not be possible to uncouple fibers 
sufficiently simply by protecting them with a cladding of sufficient 
thickness. 
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Fig. 8 — Coupling of round fibers. Only the case of horizontal polarization is 
shown. n m ka = 21, ya = 0.857, n c /n m = 1.003. 
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Fig. 9— Additional mode loss caused by the lossy surrounding medium. n m ka 
21, n c /n m = 1.003. 
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